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Abstract 

We study neutrino physics in a model with one large extra dimension. We 
assume the existence of two four-dimensional branes in the five-dimensional space- 
time, one for the ordinary particles and the other one for mirror particles, and we 
investigate neutrino masses and mixings in this scheme. Comparison of experimental 
neutrino data with the predictions of the model leads to various restrictions on the 
parameters of the model. For instance, the size of the extra dimension, R, turns 
out to be bounded from below. Cosmological considerations seem to favor a large 
R. The usual mixing schemes proposed as solutions to the solar and atmospheric 
neutrino anomalies are compatible with our model. 
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1 Introduction 



The possibility that the physical space has more than three dimensions [|I| has attracted a 
great deal of attention recently. Refs. pioneered the idea that the presence of extra 

compact spatial dimensions, some of which may be as large as fractions of millimeter, 
could lower the "fundamental" Planck scale M^, , i.e. the Planck scale in the higher- 
dimensional space-time, down to the TeV energy range, alleviating thereby the hierarchy 
problem. If the number of new dimensions is n and we assume that all of them have 
the same radii of the size R, the four-dimensional Planck scale Mpi is related to the 
fundamental Planck scale through Mpi ~ M""*"^/?". Given the experimental fact that no 
deviations from the Newton law has shown up down to i? ^ 1 mm 0], it is immediately 
clear from this relation that one extra dimension would not be enough for the solution of 
the hierarchy problem. On the other hand, for n = 2 one gets i? ~ 0.1 mm, which is in 
an interesting range in view of the gravitational force experiments discussed in 0, . 

The basic assumption of the extra dimension scenario is that all the fields charged 
under the Standard Model (SM) gauge group are localized on a brane, the familiar 3+1- 
dimensional space-time, embedded in the 4+n-dimensional space, called the bulk 
Gravitons and other particles with no SM interactions are not confined to the brane but 
are free to propagate in the bulk as well. This brane-bulk structure is suggested by some 
string theoretical observations. 

The large extra dimension scenario is particularly intriguing from the neutrino physics 
point of view. The right-handed chiral components of neutral leptons (z^/j), if they exist, 
are the only SM particles that can live in the bulk [0, ||. The Yukawa coupling between 
ul, residing on the brane, and uji is suppressed by a factor M^/Mpi^ providing a new and 
elegant explanation for the lightness of neutrinos. Perhaps the most natural conventional 
explanation is due to the seesaw mechanism, which requires the existence of a new mass 
scale around 10^^ GeV or higher. In the extra dimension scheme such a high mass- 
scale is neither needed nor naturally appears, but small neutrino masses follow from the 
suppression of the Yukawa couplings due to the large bulk volume. 

The existence of the bulk neutrinos has another interesting consequence. The bulk 
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neutrinos uin {I = e, /z, r) appear in four dimensions as Kaluza-Klein excitations having 
Dirac mass terms with the left-handed bulk components that originate in the quantized 
internal momenta in compact extra dimensions. The existence of these new degrees of 
freedom enrich the neutrino spectrum and the neutrino oscillation patterns. Theoretical 
and phenomenological aspects of this scenario have been addressed in several recent papers 
(see e.g. [i-H]). 

The localization of the SM fields on the brane may be explained in terms of non- 
perturbative effects in string theory 0. An alternative and perhaps more intuitive ap- 
proach is achieved in the context of an effective field theory. In this method one introduces 
a five- dimensional scalar field (for n = 1) with an effective potential that has a domain 



wall type profile in the extra dimension p6|-p8[. The origin of the scalar field is usually 
left unspecified. For example the so-called radion may be a natural, though not the only, 
possibility for such a scalar field. The radion field is associated to the extra dimension 
components of the metric tensor and therefore is always present in the model. 

The discussion in this paper is based on the observation that the scalar field responsible 
for the localization of the SM fermions must form two kinks, i.e. a kink and an anti-kink 
(or any number of such kink/anti-kink pairs), in order to be continuous in a compactified 
extra dimension. That kind of profile leads to interesting consequences since in addition 
to the "usual" brane, where the SM world resides, there necessarily exists at some distinct 
point of the extra dimension another brane, a mirror brane. In the mirror brane there 
live particles whose gauge interactions are identical to those of the SM particles but with 
reversed chiralities, that is, they are mirror particles Besides gravity, the ordinary 
brane and the mirror brane can communicate with each other only via possible inert 
neutrinos, that is, the right-handed ordinary neutrinos Uui and the left-handed mirror 
neutrinos Nn that live in the bulk. 

We will study in this paper the neutrino sector of the brane-mirror brane scenario. We 
neglect for simplicity any flavour mixing and consider neutrino masses and mixing within 
a single family of neutrinos and mirror neutrinos. The neutrino sector under study thus 
consists of the fields u^, u^, and Nr. In Sec. 2 a concrete realization of our scheme is 
constructed. In its framework we derive an effective Lagrangian from which the neutrino 
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mass matrix is obtained. Unlike in many previous works, our mass matrix is finite since 
we integrate out the massive Kaluza-Klein modes. In Sec. 3 we study numerically the 
generic constraints neutrino data set upon the various fundamental parameters of the 
model, such as the size of the extra dimension, the mass of the sterile neutrino the model 
predicts and the fundamental mass scale (the vev of the Higgs doublet) of the mirror 
brane. Sec. 4 is devoted to discussion of the obtained results. Concluding remarks are 
presented in Sec. 5. 



2 Theory 

Let us now construct our model for a two-brane world, where one of the branes traps the 
SM fermions and the other one the mirror fermions. The gauge inert states, the right- 
handed neutrino uji and the left-handed mirror neutrino Nl, are assumed to propagate 
freely in the whole higher-dimensional bulk. For the localization mechanism of the gauge- 



active fermions, based on effective field theory approach, we follow Refs. For 
simplicity and clarity we restrict ourselves in the following to the case where there is 
only one extra dimension, but the extension of the analysis to the cases of several extra 
dimensions is quite straightforward. 

We assume the extra dimension to be compactified on a circle of radius R. The 
coordinate system reads z = (x^, y), where /x = 0, 1, 2, 3 and y ~ y + 2TTR. The scalar field 
responsible for the localization of fermions on the 4-dimensional membranes is denoted by 
$. It is assumed to have two kinks, one at y = where its value grows from — / to +/ and 
another one a.t y = where its value drops back to — /. Between the kinks $ is supposed 
to remain constant. Fermions are localized on a finite- width wall around the points in the 
fifth dimension where $ = as, heuristically speaking, their position-dependent masses 
are there the smallest. Outside these positions field fluctuations are strongly suppressed 



due to higher masses [28 



The five- dimensional spinors are decomposed as 



L 
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where z/^ and z/^ are related to the ordinary neutrinos and A^^ and N'j^ to mirror neutrinos. 
These fields are functions of the five-dimensional space-time. The appropriateness of our 
notation will become clear towards the end of this section. 

We assume the relevant part of the brane-world action to be of the form 

S = J d^xdy[L{iV^dA -m + ^)L + i^T^S^^ + {kH*^L + h.c.)] , (2) 

where A = 0, . . . , 4, m is a mass parameter, k is a dimensionful Yukawa coupling, and H 
is the neutral component of the standard SU(2)-doublet Higgs field. The gamma matrices 
are given by 

and all the fields are taken to be periodically continuous, e.g. L(a;, 0) = L{x,2nR). The 
term L$L takes care of the localization of the active neutrinos on the branes. In general 
the action could also include Dirac and Majorana mass terms for bulk neutrinos but 



we neglect them in the present discussion. It should also be remarked that the presence 
of the inert neutrinos in the bulk is an assumption, because we have omitted in the action 
the term 'l'$\I', which would have localized the inert neutrinos on the branes. 
The equation of motion for L is 

iV^dAL = {m- $)L , (3) 

where a term k*H^ has been neglected as in classical approximation \E' = 0. Supposing 
that $ = $(?/), and writing A''^ = Nji{x)gR{y) and u'j^ = i^Lix)gL{y), one ends up with 

igL{y)(y^'d^uL = Nn^dy + m- ^y))gR{y) , 
i9Riy)cT''d^NR = ULi-dy + m- ^y))gLiy) ■ (4) 

The left-hand side of Eqs. (Q) being classically zero, it is trivial to see that gL,R{y) have 
solutions behaving so that gR^y) tends to confine near y = y^ and giiy) near ?/ = 0. 
In other words, since and Nr really get localized on different branes, the chosen 
action reproduces exactly the features that we have been looking for. The consistency 
of the theory requires that Jq^^ dy{^{y) — m) = 0, which for thin branes reduces to 



5 



71 — Ad = 7rm//, where A6 = 27i — y^/R. Slightly more contrived considerations, which 
we shall not present here, reveal that the Higgs field has a twin-peaked profile with 
maxima at y = and y = y^. One may thus approximate (somewhat symbolically, cf. the 
appendix of Ref. |]^) 

y'L = ^)yL{x), N'^ = ^5{y-y.)NR{x), H = ^)h.{x) + ^5{y - y.)h+{x) , (5) 

where h-{h^) can be viewed as a vev of the Higgs field in "our" brane (the mirror brane). 
By substituting these, together with the familiar Kaluza-Klein expansion 



V N'L{x,y) J 
to the original action, one has finally 



E 



Jny/R 



(6) 



I n=— oo 



(7) 



where 



u 



and K is taken to be real. 

The first six ra-dependent terms of the action (|^) disappear due to the equations of 
motion of the Kaluza-Klein excitations z/^j^ and Nin (with n 0), 

iR 



T^Rn 



n 



iR 



n 



ia^d.NLn + uhlNRe 



inAe\ 



(9) 



The remaining two terms of the action then yield 



E 



n J n 



R, 



R 



-u—h+N\^a^'d^VRne ' 



inAe 



n 



(10) 
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or, by applying Eqs. (||) again, 

£„^o = u^siAe)Rih^h*_NliyL + h.c.) + Y.iu''^h_ula''d^ih*_UL) 

+ E tu'^h^Nla'-d.ihlNR) , (11) 
where terms including higher derivatives of v^n and Nin have been neglected, and 

s{Ae) = i = TT - (A^ ^ 0) . (12) 

This on-shell Lagrangian could be equivalently obtained by integrating out the massive 
Kaluza-Klein excitations. 

Note that since h\. are constants in the leading order, they can well be taken out of the 
derivatives in Eq. ([TT]). Examinations show that the values of h± are essentially dependent 
on the widths of the ordinary and mirror branes. The quantity /i+ can be considered 
here as a free parameter while h- is bound by the usual Higgs scalar expectation value 
= 174 GeV. 

Combining finally Eq. (pAj) and Eq. (^ (for n = 0) with suitable rescalings of the fields 
so that the kinetic terms take the canonical form, one ends up with a Lagrangian 

+u'^s{Ae)R{n+n_h+h*_NliyL + h.c.) + u{n_h*_v^j^QVL + n+hlNl^Nn + h.c.) , 

where 

1 

and the identity Z^n^o^^^ = vr^/S has been used. In this Lagrangian, where all the 
massive Kaluza-Klein modes are integrated out, only the lowest order effects of the extra 
dimension {i.e. field rescalings and mass terms) are kept. 
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3 Neutrino phenomenology 



According to the Lagrangian L^tt the neutrino mass term of the model is given by 



M 



-{VlVrNrNl) 
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(14) 



where 



un\h^ 



b = u'^s{A9)Rn+n^\h+h. 



c = un^ I Ha 



(15) 



with appropriate redefinitions of the fields, and we have associated vr and with the 
zero modes of the respective Kaluza-Klein states. More familiarly the mass Lagrangian 
can be presented in the form 



M 



-^rMMl + h.c. 



(16) 



where N"[ = (z/^ z/£ N^), and M is the matrix appearing in Eq. ([T^). 

As one can see from the results above, the dependence of the neutrino mass Lagrangian 
on the fundamental parameters of the theory, such as the radius R, the values of the scalar, 
h±, the coupling k, and the relative positions of the brane and mirror brane {i.e. A6), is 
quite non-trivial and difficult to analyse analytically. We will therefore study numerically 
the generic constraints the present neutrino data sets on the model by varying unknown 
parameters within conceivable range of values and plotting the ensuing predictions for 
various measurable quantities. 

Let us start the phenomenological analysis by diagonalizing the neutrino mass matrix 
obtained above. Since M is real and symmetric, the diagonalization can be performed 
by an orthogonal matrix O. (We consider here just one neutrino family. The different 
families may have different mixing angles.) 

If we define 

/ sin a sin a cos a — cos a \ 



O 



1 

V2 



cos a 

— sin (3 sin /5 — cos /5 — cos /5 

cos /3 cos (3 sin /3 sin /3 

V cos a cos a — sin a sin a ) 



(17) 
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with 



sin^ a 



we obtain 



+ 



+ 6^ - 



2 2 J(a2 + 62 + c2)2 - 4a2c2 



sm^ /? = - + ^ 



(18) 



O^MO 



( -m+ \ / ^1"^ 

m+ 

-m_ 

V m- J 



\ 

(T2m+ 

(T3m„ 

V a4m_ / 



(19) 



where the eigenvalues are 



a^ + b'^ + c'^± Jia^ + 62 + c2)2 - Aa^c^ 



(20) 



and (Tj's are sign factors. Note that if 6 = 0, matrix O can be presented in a much simpler 
form. This case, however, is realized only in the special case of AO = it, i.e. when the 
two branes are in opposite locations in the fifth dimension. 
The mass eigenstates are given by 



in terms of which the mass Lagrangian reads 



(21) 



i™.(xa. + x.x.)-i™-fex-3 + x-*). 



By defining 
one has finally 



V' = ^(Xl+X2), 0=-^(X3-X4) 



:22) 



(23) 



(24) 



Cm = —m+ipip — m_(j)(j) . 

The neutrino sector thus consists of two Dirac neutrinos, as was expected as the theory 
conserves, by construction, lepton number. 

From Eqs. ([T7|) , ( pT]) and ( p3D it follows that the ordinary left-handed neutrino vi is 
the following superposition of the mass eigenstate neutrinos: 



vl = cos a (pL + sin a ip^. 



The superposition orthogonal to this combination is the inert field N^. As ul is the only 
active neutrino living in our brane, only the mixing angle a can be experimentally probed. 
The other angle (3 parametrizes the mixing between the right-handed fields and N^i, 
and is therefore a measurable quantity only in the mirror world. 

There exist several empirical constraints on the masses and mixing angles of neutri- 
nos, coming from laboratory experiments, astrophysical observations and cosmological 
considerations. In our case, where fiavour mixings are neglected, mixings occur between 
an active and a sterile neutrino of each family. It turns out that the existing laboratory 
bounds on the mixing angle a between the electron neutrino and a heavier neutrino in 
various mass ranges of the heavier neutrino mass are in general ineffective in our model. 
This is because the mass of the extra neutrino turns out to be quite low for a plausible 
choice of the model parameters. The only relevant laboratory constraints are due to the 



electron neutrino mass measurement in the tritium beta decay experiments |30| and the 
neutrino oscillation experiments. The upper limit of the electron neutrino mass from the 
beta decay experiment is [^] 

m,, < 2.3 eV. (25) 

As far as the electron neutrino is concerned the most stringent oscillation limit comes from 
the Bugey disappearance experiment PT|. This limit has been approximated as follows 
in our numerical calculation: 

sin22a<0.1 for 100 eV^ > l^m^j > 2 eV^ , 

sin2 2a<0.02 for 2 eV^ > l^m^j > 0.04 eV^ , (26) 
sin2 2a < 0.1 for 0.04 eV^ > \5m'^\ > 0.01 eV^ . 

The active-sterile mixing can also be constrained by cosmological arguments. If the 

mixing is too large, neutrino oscillations, acting as an effective interaction, would bring 

the sterile neutrino in equilibrium before neutrino decoupling, and the resulting excess 

in energy density would endanger the standard scheme for the nucleosynthesis of light 

elements (BBN) [0. This leads to the following bound for ^ Vs mixing [Q: 

\5m?\ sin^ 2a < 5 x 10^^ eV^ for \5m'^\ < 4 eV^ , , . 

sin^ 2a < 10"* for j^m^j > 4 eV^ . ^ ^ 

This bound is, however, avoided if there was a suitable net lepton number in the early 



universe 34 . 
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With applying these constraints, we search the allowed regions of the parameter space 
numerically using a Monte Carlo analysis where we have varied three unknown parameters: 
the radius of the extra dimension i?, the mirror brane Higgs value /i+ and the mirror brane 
position parameter The extra dimension radius is varied from the Planck scale to 
a millimeter scale and the mirror brane Higgs value between 10^^ MeV and 10^° MeV. 
We consider these parameter ranges wide enough for a representative analysis. Both R 
and are randomized so that their logarithms are evenly distributed. For the position 
parameter A6 an even distribution between and vr has been taken. 

In addition to these parameters the action depends on the dimensionful coupling k. 
The natural scale for it could possibly be set by the higher- dimensional Planck scale M^,, 
and therefore we write (with Mpi = M^R) 

1/6 



K = ^= = k' , (28) 

^/M: iMi ' ' ^ ' 



k' , I R 

where k' is a dimensionless coupling constant. Its value is not really known, but a plausible 
choice would be a number relatively close to unity. We have performed our analysis for 
the values k' = 1 and = 0.01. 

A central parameter in the scheme is obviously the size R of the extra dimension. 
From the experiments that probe the effects of gravity in short distances one knows that 
i 

R < 1 mm. (29) 

From the dependence of the elements of the mass matrix M on R follows an upper bound 
on the larger mass eigenvalue m_|_ as a function of R, as presented in Fig. 1. When the 
limit ( |2^ ) is saturated, the largest allowed value is m_|_ ~ 10 eV, but larger values are 
allowed for a smaller R. 

4 Results and discussion 

Let us now proceed to our numerical results. In Fig. 2 we present a scatter plot in 
the [R, sin a;)-plane, obtained by allowing the model parameters vary as indicated above, 
with and without taking the cosmological constraint, Eq. (pT]), into account (marked in 
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logfi^m) 

mm 



Figure 1: Maximal neutrino mass m™*^^ (dashed line) and coupling k! (solid line) 
as function of minimal radius of the extra dimension i?min- 



the plot by crosses and dots, respectively). Fig. 2a corresponds to the case k' = 0.01. 
Let us now look at the general features displayed by this figure. Comparison of Eqs. (H), 



(p!5|) and ( pSf ) implies that a,c ^ b for small values of R. From Eq. (|T8]) it can then be 
inferred that, depending on the relative sizes of and h-, either sin^ a ~ 1 (/i+ < /i_) 
or sin^ a ~ (/i+ > h^). In other words, for small values of the size R of the extra 
dimension, the mixing angle a has values either close to or close to 7r/2, the other 
values being forbidden. The small values of a correspond to the case where a < c, that 
is, where the ordinary active neutrino vl is predominantly the lighter mass eigenstate 0^ 
and the sterile mirror neutrino A^^^ is predominantly the heavier state ipi. The values of 
a close to ti/2 in turn correspond to the situation where a > c, and vl is predominantly 
the heavier state ipi and the sterile mirror neutrino A^^^ is predominantly the lighter state 
0L, the both neutrino masses now being below the experimental limit m± ^ 2.3 eV. 

It is important to note that the size of the extra dimension, R, is bounded from below 
in our scenario. This is an obvious consequence of the fact that the quantities a and c 
increase with decreasing i?, and the mass of the predominantly active neutrino state ~ vl-. 
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log(-^) 

mm ' 



Figure 2: The neutrino mixing angle sin (a) as function of extra dimension radius 
R for k' = 0.01 (a) and k' = 1.0 (b). Crosses indicete points allowed both oscillations 
and cosmology, whereas dot are points allowed by oscillations only. 

whose mass experiments test, is determined mainly by these quantities, since for small R 
a,c^ b. In the case of k' = 0.01(1.0) the lower limit is i? ^ 4 x 10"''(5 x 10~^) mm. 
More generally the lower limit of R is obtained from 

^^-^ < {RMpif/^ X 2.3 eV , (30) 
V 27r 

as long as k' ^ 0.1 (c/. the sohd line in Fig. 1). 

When the cosmological constraint, Eq. (|27|) , is taken into account (crosses in Figs. 2a 
and 2b), large values of R are favoured (-R ^ 3 x 10^^ mm), unless the mixing angle a is 
very close to or 7r/2. This is because for smaller R the squared mass difference — 
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tends to increase, jeopardizing the fulfillment of the cosmological condition. Let us note 
that in the parameter space covered by our Monte Carlo analysis, the cases allowed by 
the cosmological condition are concentrated to the region of large mixing angle a, i.e. to 
the case where the active neutrino is heavier than the sterile one. Nevertheless, also the 
small values of the mixing angle are allowed, i.e. the case where the sterile neutrino can 
be relatively heavy, albeit with quite a small part of the parameter space. 

The step-like behavior of the scatter plot in the region 10~^ mm ^ R ^ 10~^ mm 
is due to the oscillation constraint, Eq. (^61). As remarked before, this bound together 
with the upper bound of 2.3 eV of the electron neutrino mass, is the only laboratory 
constraint that is effective for our model. The reason is the relatively small values of the 
inert neutrino mass the theory allows, as displayed by Fig. 1. 

Fig. 2b is the same scatter plot as in Fig. 2a, except that now n' = 1. A comparison 
with the previous case reveals the interesting fact, namely that for a large values of the 
coupling k' only the large mixing angles a are allowed. This can be understood as follows. 
The ratio of the mass matrix elements a and c, 



a n_\h-\ \h_ 



l + 7r%2i?2|/^^|2/3 ^^^^ 



c n+\h+\ 1 + 7r%2i?2|/^_|2/3 ' 

is approximately equal to 1 when the product u^R^ is large enough. From Eq. ([T8| ) one 
can deduce that in this case sin^ a ^1/2. In the case k' = 1 one is practically always 
in this regime, taken into account the lower limit for R from the electron neutrino mass 
measurement. In the previous case of n' = 0.01 this regime is reached, because of the 
smaller value of u, only at the upper end of the allowed R range, as can been seen in Fig. 
2a. The lower bound for R, still given by the experimental upper limit for the electron 
mass, is now larger than in the case of k' = 0.01 since the neutrino mass (for small b) 
is proportional to k', as seen from Eqs. (|T3p, ( PP] ) and When R is increased, all 

possibilities from ul — ipL to the maximal mixing ul = {4>l + '^l)/^/^ become open. 
At low R cosmology forces vl to be very close to the heavier mass state ipi, but when 
i? ^ 3 X 10~2 mm, may develop a large component as well. 

Fig. 3 presents the region in the parameter space (sin^ 2a, 5rn?) (with 5rn? = m\ — im?_) 
reached by our model in the case k' = 1. The only constraint taken into account here is 
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Figure 3: Neutrino mass difference as a function of sin^(2a) . The solid line divides 
the space to excluded (upper part) and allowed (lower part) regions by cosmological 
considerations {k' = 1.0). 

the upper bound 2.3 eV for electron neutrino mass. We have checked that the inclusion of 
other laboratory bounds do alter Fig. 3, but only cuts ~ 1% more points from the region 
which is anyway excluded by oscillation experiments. The line between the cosmologically 
allowed (lower part) and excluded (upper part) regions has been drawn to Fig. 3, too. 
For k' — 0.01 the allowed region would extend to larger mass differences, with virtually 
no change in cosmologically allowed region. As one can see, our model allows for all 
the active-sterile mixing schemes discussed as possible solutions to the solar neutrino 
and atmospheric neutrino anomahes. In particular the so-called SMA solution (|5m^| 
(3 - 10) X 10"^ eV^, sin^2Q; fa (0.2 - 1.3) x 10~^) of the solar neutrino anomaly can 
be realized in this model. Note, however, that the solar LMA solution is excluded by 
cosmology bound. 
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5 Conclusions 



We have investigated a five-dimensional model with two four-dimensional branes, a brane 
for ordinary SM particles and another brane for mirror particles. By looking at neutrino 
masses and mixings, we have shown that confrontation of the predictions of the model 
with data constrains the extra dimension physics remarkably. The experimental upper 
bounds on the masses of the known neutrinos directly restrict the extra dimension size, 
independently of how neutrinos actually mix. The value of brane-bulk neutrino coupling k' 
determines the minimal value of extra dimension size i?, a larger coupling corresponding 
to a larger -Rmin- The nature of neutrino mixing depends crucially on the coupling n' . 
For k' ^ 1 the predominantly active neutrino, i.e. the ordinary neutrino, is never lighter 
than the predominantly sterile neutrino. For smaller values, k' ~ 0.01, also the opposite 
situation is allowed. In both cases active and sterile neutrino can mix maximally when R 
is large enough. 

As may be seen from Fig. 1, also the largest possible mass of the heavier neutrino 
(of each flavour) depends essentially on k' , a smaller k' {i.e. a smaller -Rmin) enabling m™^ 
to be higher. For k' = 0.01 the maximal neutrino mass is ~ 1 keV. However, already 
for k' ~ 10~^ the neutrino mass may be of the order of 1 MeV with -Rmin ~ 10'^*^ mm. 
Therefore the issue of the value of k' is emphasized. 

If also cosmological constraints are taken into account, the room for allowed neutrino 
mixings is largely suppressed. With cosmological bound the model clearly favor large 
radius -R ^ 10~^ mm. Only mixing angles very close to 7r/2 (and equally to for some 
values of k') make exception, all the radii down to the minimal one being then allowed. 
In this case no direct method exists for observing mirror brane neutrinos. 

The model we have studied is the simplest possible in the sense that only one extra 
dimension is introduced. This kind of approach is not fully realistic since it is known that 
the hierarchy problem cannot be accounted for unless the number of extra dimensions is 
at least two. It is clear that with more than one extra dimension the brane-bulk structure 
becomes much more complicated than in the five-dimensional case we have studied in this 
paper, increasing the flexibility in the values of the model parameters. Nevertheless, we 
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expect that our model exhibits the generic features that also appear in higher-dimensional 
treatments. 
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